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I.  INTRODUCTION 


Shaped  charge  warheads  are  utilized  in  a large  percentage  of 
military  amiuanitions . The  processes  of  warhead  design,  development 
and  testing  are  rather  lengthy,  expensive  and  ^ime  consuming. 
Mathematical  modelling  and  computer  simulation  of  the  performance  of 
warhead  designs  can  reduce  considerably  the  amount  of  proof  testing 
required  and  eliminate  unnecessary  development  costs.  Parametric 
studies  of  certain  configurations  or  of  proposed  design  changes  can 
be  performed  quickly  and  economically  once  a basic  model  has  been 
validated. 

Most  conventional  shaped  charge  warheads  consist  of  a metallic 
conical  liner  with  a vertex  angle  of  30“  to  90“ , and  a cylindrical 
charge  which  is  molded  around  it  and  is  point  initiated  at  its  end. 

The  resulting  detonation  wave  is  approximately  planar  by  the  time  it 
strikes  the  liner  and  causes  the  liner  to  collapse  resulting  in  a low 
mass  high  speed  jet  moving  forward  and  large  mass  low  speed  slug 
moving  behind  the  jet.  The  mass  partition  of  the  different  liner 
elements  into  jet  and  slug  portions  takes  place  at  their  respective 
stagnation  points.  The  theory  of  jet  formation  for  such  warheads  was 
first  published  in  the  open  literature  by  Birkhoff,  MacDougall,  Pugh 
and  Taylor^  and  was  later  extended  by  Pugh,  Eichelberger  and  Rostoker^. 
The  whole  theory  was  based  on  simulating  the  relative  motion  of  the 
metallic  liner  during  collapse  with  the  steady  flow  of  two  jets  of 
water  impinging  upon  each  other  at  the  stagnation  point.  The  applica- 
tion of  the  above  theory  necessitated  assuming  the  collapse  velocity 
of  each  liner  element  or  measuring  it.  Kiwan  and  Wisnieski  enhanced 
the  above  theory^  by  calculating  numerically  the  collapse  velocities, 
jet  and  slug  characteristics  from  the  explosive  properties  and  the 
assumed  geometry.  They  demonstrated  the  procedure  by  calculating  the 
properties  of  two  wedge  shaped  charge  liners. 

Recent  interest  has  been  shown  in  exploiting  other  technologies 
to  develop  new  shaped  charge  warhead  designs  for  future  weapons 
of  the  1980s.  In  this  report  we  shall  investigate  some  designs  of 
interest  and  study  their  characteristics  through  numerical  simulation 
of  their  performance.  The  design  shown  in  Figure  1 consists  of 
a hemispherical  metallic  copper  liner  and  a hemispherical  charge 


^G,  Birkhoff j D.  P,  MaoDougallj  E.  M,  Pugh,  and  Sir  G,  I.  Taylor 
"Explosives  with  Lined  Cavities"  J.  Appl.  Phys.  19(2948)  p.  563. 

^E.  M.  Pugh,  R,  J.  Eichelberger,  and  N.  Rostoker  "Theory  of  Jet  For- 
mation by  Charges  with  Lined  Conioal  Cavities"  J.  Appl  Phys.  23, 
(1961)  p.  532. 

^A,  R,  Kiwan  and  U.  Wisniewski,  "Theory  and  Computations  of  Collapse 
and  Jet  Velocities  of  Metallic  Shaped  Charge  Liners"  BRL  Report  No. 
1620. 
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together  with  an  initiation  package.  The  major  portion  of  this  study 
will  deal  with  tho  hydrodynamic  simulation  of  liner  collapse,  jet 
formation,  flight  and  characteristics  obtained  from  the  charge  shown 
in  Figure  1.  Figure  2 shows  an  idealized  schematic  of  the  anticipated 
liner  collapse  after  being  hit  by  a convergent  detonation  wave.  A 
typical  element  b collapses  towards  the  center  C with  velocity 
after  being  hit  by  the  convergent  detonation  wave.  An  observer  at 
the  pole  P sees  the  element  E collapsing  towards  him  with  velocity  V^. 
The  liner  continually  get  thicker  while  it  is  being  compressed,  until 
eventually  it  starts  jetting.  The  dotted  contour  shows  the  actual 
liner  configuration  at  that  time.  The  equatorial  section  of  the 
liner  is  elongating  due  to  the  expansion  of  the  detonation  products  at 
the  unconfined  equator. 


II.  HYUROUYNAMIC  CONSIDERATIONS 

The  hydrodynamic  simulation  of  the  collapse,  jet  formation  and 
characteristics  of  the  above  referenced  charge  is  made  computationally 
on  the  generic  charge  shown  in  Figure  3 which  does  not  contain  the 
initiation  package  referenced  in  Figure  1.  Only  half  of  the  charge 
is  shown  in  Figure  3 due  to  it  being  axisymmetric.  The  HELP  code** 
was  employed  to  simulate  computationally  the  performance  of  the  charge 
shown  in  Figure  3. 

HELP  is  a two  dimensional  finite  difference  multi  -material 
Eulerian  code  capable  of  treating  compressible  fluids  and  solids  in 
the  hydro  and  elastic  plastic  regimes.  The  conservation  equations 
that  are  solved  in  HELP  are: 

3p  

at  " “ ax^ 

^ Dt  ' “ ax^ 

^^T  = a 

" ^ 3x7 

where  Xi  denotes  the  i^^  Goordi.iar-*  of  position,  t the  time,  p denotes 
the  density,  ui  the  i^^  veloci.y  .proponent,  E-j-  the  total  energy. 

D stands  for  the  total  material  derivative,  a.,  is  the  stress  tensor 
Dt 


^L,  J.  Uagman  et  at,,  "HELP  A Multi  -material  Eulerian  Program  for 
Compressible  Fluid  and  Elastic  - Plastic  Flows  in  Two  Space 
Dimensions  and  Time"  Systemst  Science  and  Software  Report  75-2654  (1975), 
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AXIS  OF  SYMMETRY  (CM) 


which  can  be  separated  into  its  hydrostatic  stress  component  (pressure) 
“ and  a stress  deviator  tensor  Sij»  i.e., 


'ij 


hj  - hj  P 


C4) 


Equations  (1) , (2)  and  (3)  are  transformed  into  conservation  form^ 
and  then  put  into  an  equivalent  integral  form  before  applying  finite 
difference  methods  to  these  equations.  The  equations  that  are  therefore 
solved  by  finite  difference  methods  for  each  cell  of  the  computational 
mesh  and  every  computational  cycle  are: 


Am 


ACm  Uj)  = 


ACmE^) 


- At 

S 

/S..n.dS  - At/pu.dS  - At/pu.u.n.dS 
iji  J 3 1 3 T- 

S S S 

/S.  u.n.dS  - At  fpu.n.dS  - At/pu.E_n.' 

iJJi  J ^ ^ J - ^ ^ 


(5) 


(6) 


,dS.  (7) 


n.  denotes  the  tmit  vector  normal  to  a surface  element  dS. 

1 

Equations  (5)  through  (7)  are  supplemented  by  appropriate  initial  and 
boundary  conditions  to  initiate  the  flow.  Material  properties  are 
represented  in  HELP  by  the  Tillotson  Equation^  of  state  for  inert 
materials,  the  ideal  gas  equation  of  state,  and  the  JWL  equation  of 
state  for  explosion  products.  The  code  also  contains  a bum  routine 
based  on  the  JWL  Equation  of  state.  The  above  flow  equations  are 
integrated  in  three  phases  corresponding  to  the  evaluation  of  the 
three  different  types  of  integrals  occuring  in  these  equations  and 
called  the  SPHASE  (Strength  phase),  the  HPHASE  (Hydro  phase),  and 
the  TPHASE  (Transport  phase). 

HELP  contains  a variety  of  options  such  as  transmittive  and 
reflective  boundary  conditions  at  various  grid  boundaries.  It  contains 
also  an  artificial  viscosity  option.  Material  interfaces  are  defined 
in  Lagrangian  manner  and  are  identified  by  massless  tracer  particles. 

A slide  line  can  also  be  introduced  along  a material  interface  separa- 
ting two  materials  which  permit  them  to  slide  against  each  other. 


®P.  D.  Lax,  "Weak  Solutions  of  nonlinear  Hyperbolio  Equations  and 
Their  Numeriaal  Computation",  Comm,  on  Eure  and  Applied  Math, 

Vol  III,  p 159  - 193,  (1954). 

SJ.  H.  Tillotson,  "Metallia  Equations  of  State  for  Hypervelooity 
Impact,"  General  Atomic  Report  GA~S216  (1962). 


11 


Although  HKLP  models  the  strength  properties  of  materials  and  contains 
n yield  criterion  and  a failure  criterion  our  calculations  will  be 
restricted  to  the  hydrodynamic  phase,  with  only  viscous  effects 
incorporated.  Our  experience  with  calculations  incorporating  the  SPIIASB 
of  the  code  is  limited,  and  we  have  not  investigated  the  validity  of  the 
model  used  in  that  part  of  the  HhLP  code.  Walsh‘S  reported  that  the 
results  of  his  calculations  using  the  HhLP  code  in  similar  problems  were 
not  significantly  affected  when  the  strength  effects  wore  incorporated. 


III.  COMl’UTATIONAL  SET  UP  AND  SIMULATION 

The  computational  simulation  of  the  performance  of  the  charge  shown 
in  Figure  3 was  made  in  a computational  mesh  of  50  x 90  cells.  Each 
cell  was  0,8  mm  x 0.8  mm  in  the  region  containing  the  metallic  liner 
whose  thickness  is  1.9  mm.  The  zones  were  gradually  enlarged  radially 
and  axially  to  the  end  of  the  mesh.  The  choice  of  the  mesh  size  had  to 
be  balanced  between  the  desire  for  accuracy  of  the  solution  obtained 
and  the  speed  of  the  computations.  A coarse  mesh  size  is  detrimental 
to  the  accuracy  of  the  calculations  particularly  so  if  any  material 
region  consisted  of  mixed  cells  only.  A fine  mosh  enhances  the  com- 
putational accuracy  but  requires  a larger  mesh  and  reduces  the  magni- 
tude of  the  time  step  and  thus  requires  longer  computational  time  to 
solve  a given  problem.  Ten  equally  spaced  latitude  circles  (initiation 
rings)  were  selected  on  the  outer  hemispherical  surface  of  the  charge 
shown  in  Figure  3.  The  simultaneous  initiation  of  these  rings  was 
considered  to  simulate  reasonably  well  the  initiation  process  of  the 
actual  charge. 

Figure  4 shows  the  liner  configuration  at  t=].92  ps  together  with 
the  velocity  field.  The  liner  is  seen  starting  to  collapse  after 
being  hit  by  the  almost  convergent  detonation  wave  about  0.46  ps  earlier. 
The  lack  of  confinement  on  the  equatorial  plane  allows  the  detonation 
products  to  expand  rapidly  from  that  surface  causing  a departure  from 
the  idealized  collapse  depicted  in  Figure  2.  The  equatorial  section  of 
the  liner  starts  to  elongate  at  the  explosive-metal-air-interface. 

As  the  equatorial  rarefaction  wave  travels  along  the  surface  of  tile 
liner  towards  the  polar  region  the  pressures  and  velocities  are  reduced. 
Figure  5 shows  the  pressure  in  the  flow  field  which  is  found  to  be 
large  and  rising  in  the  collapsing  liner  where  it  has  a maximum  value 
of  0.328  Mbar.  About  that  time  (t=1.92  ps)  the  effect  of  the 
rarefaction  from  the  nearest  unconfined  spherical  surface  of  the  charge 
starts  to  be  felt  at  the  liner  surface  a.nd  the  pressure  begins  to 
decrease.  Figure  6 shows  the  flow  field  at  t^s.S  ps.  The  liner  is 
observed  to  be  getting  thicker  at  the  pole  and  elongating  further  in 
the  equatorial  region.  The  pressure  field  is  seen  in  Figure  7 to  have 


V.  M.  Walsht  private  oormunioation. 
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Figure  5.  Pressure  Field  at  an  Early  Stage 
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decreased  significantly  by  that  time.  The  maximum  pressure  in  the 
liner  material  is  found  to  be  O.IS  Mbar.  As  the  liner  collapse 
advances,  the  pressure  starts  to  increase  again  due  to  liner  com- 
pression and  reaches  a maximum  value  of  0.603  K&ar,  at  t=9.2  us 
as  can  be  seen  in  Figure  8,  The  pressure  decreases  thereafter  due 
to  the  influence  of  the  equatorial  rarefaction  and  the  expansion  of 
liner  material  arising  from  jet  formation  and  elongation.  The 
rarefaction  wave  emanating  from  the  equatorial  plane  reaches  the  polar 
region  about  t-7.2  us  and  the  jet  becomes  distinguishable  after  that 
time.  Figures  9 and  10  show  the  early  stages  of  jet  formation.  The 
short  arrows  show  the  direction  of  the  local  flow  velocity  in  the  differ- 
ent layers  of  liner  material.  Figures  11  through  13  show  the  late 
stages  of  jet  formation  and  flight.  A particularly  remarkable  feature 
in  those  Figures  is  that  the  equatorial  section  of  the  liner  impacts 
the  jet  after  its  initial  elongation.  Several  other  features  of  these 
jets  will  be  discussed  more  fully  later.  Figure  14  shows  a plot  of  the 
velocity  versus  time  for  ten  different  tracer  particles  placed  on  the 
inside  surface  of  the  hemispherical  liner,  while  figures  15  through  17 
show  plots  of  velocities  versus  time  of  four  particles  placed  across  the 
liner  thickness  at  throe  different  locations,  one  on  the  axis  of  symmetry 
and  the  others  are  placed  at  locations  whose  radii  make  angles  of  36“ 
and  72“  respectively  with  the  axis  of  symmetry. 


IV.  RESULTS  AND  EXPERIMENTAL  COMPARISON 

The  primary  advantage  of  mathematical  modelling  and  computational 
simulation  of  physical  problems  is  the  wealth  of  information  available 
6nce  a successful  model  has  been  achieved.  The  values  of  the  various 
calculated  physical  parameters  of  the  problem  are  available  throughout 
the  region  of  computation  in  a permanent  record  form  which  can  subse- 
quently be  retrieved  and  examined.  Parametric  studies  can  also  be  made 
quickly  and  sometimes  economically.  In  this  section  we  shall  summarize 
our  computational  results  and  compare  them  with  some  of  the  available 
experimental  data.  Figure  18  shows  a plot  of  the  average  velocity 
components  of  the  metallic  liner  as  a function  of  time,  while  Figure  19 
shows  the  different  liner  energies  as  functions  of  time.  The  initial 
rarefaction  wave  arriving  at  the  liner  surface  due  to  the  nonconfinement 
of  the  spherical  charge  surface  reduces  the  liner  acceleration,  while 
the  equatorial  rarefaction  wave  causes  the  jet  to  become  distinguishable. 
The  continued  liner  compression  converts  the  liner  radial  momentum  to 
axial  momentum.  The  total  liner  energy  increases  rapidly  at  first  and 
approaches  an  asymptotic  value  later  on.  Figure  20  shows  a plot  of 
jet  velocity  as  a function  of  cumulative  mass  at  various  times.  The 
jet  mass  continually  increases  as  more  metal  is  accelerated  to  jet 
velocities. 
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It  was  found  numerically  that  about  18.7%  of  the  hemispherical 
liner  forms  the  jet  (i.e.,  has  velocity  ^2nim/ijs).  lixperimental 
measurements  estimate  the  jet  mass  to  be  about  18%  of  the  liner  mass. 

The  calculations  indicate  that  about  19.9%  of  the  explosive  energy 
is  delivered  to  the  copper  liner  of  which  11.8%  is  in  kinetic  energy 
form  and  8.1%  is  in  internal  energy  form.  The  kinetic  energy  of 
the  jet  is  41%  of  the  total  energy  of  the  liner  (69%  of  the  kinetic 
energy  of  the  liner).  The  jet  tip  velocity  was  found  to  be  6.42mm/ps. 

Two  experimental  measurements  of  jet  tip  velocity  wore  made  at  BRL  of 
7.07mm/us  and  7.57mm/us.  In  the  experimontal  tests  the  hemispherical 
charge  used  was  0.151  kg  of  PBX  and  the  initiation  package  contained 
0.2685  kg  of  composition  B-3.  In  the  computations  only  the  energy 
from  the  hemispherical  charge  was  incorporated  in  the  calculations . 

If  one  adds  the  kinetic  energy  of  the  flier  plate,  in  the  initia- 
tion package,  to  the  explosive  energy  of  the  charge  then  the  resulting 
jet  velocities  will  increase  on  the  average  by  about  8%.  The  resultant 
jet  tip  velocity  will  increase  to  6.92mm/vis.  If  in  addition  one 
considers  the  effect  of  the  confinement  provided  by  the  flier  plate 
and  scale  the  jet  tip  velocity  according  to  the  values  given  in  columns 
(a)  and  (e)  of  Table  1 of  this  report,  the  jet  tip  velocity  will  be 
found  to  be  7.35mm/us  which  is  within  the  range  of  experimental  measure- 
ments . 

Several  important  observations  can  be  made  from  the  computations , 
concerning  the  jet  properties  from  such  systems.  No  inverse  velocity 
gradient  was  found  in  the  jet  produced  by  such  systems  as  in  the  case 
of  most  conical  systems.  The  lead  particle  in  the  jet  has  as  a result 
of  this  a small  amount  of  mass.  Examination  of  Figures  9 and  10  suggests 
that  the  jet  forms  from  the  liner  material  located  on  the  inside 
surface  layer  of  the  liner,  while  the  remaining  liner  material  ends 
up  in  the  slug.  This  observation  is  confirmed  from  an  examination  of 
the  velocity  plots  shown  in  Figures  14  through  17.  Figure  14  shows 
that  the  entire  inside  surface  layer  of  the  liner  reaches  terminal  jet 
velocities  except  for  the  equatorial  portion  of  that  layer.  Figures 
15  through  17  confirm  this  observation.  The  amount  of  jet  mass  and 
kinetic  energy  from  such  a system  is  similar  to  that  obtained  from 
conical  systems  (same  order  of  magnitude) , although  their  distributions 
along  the  jet  length  appear  to  be  different.  X-rays  from  early 
experimental  tests  showed  a part  of  the  jet  to  be  missing.  It  was 
thought  at  the  time  that  that  part  of  the  copper  jet  vaporized.  The 
calculations  shown  in  Figures  11  and  12  show  the  correct  interpretation. 
It  is  clear  from  those  figures  that  the  equatorial  portion  of  the 
hemispherical  liner  after  initially  stretching,  moves  inwards  and 
pinches  off  the  jet.  X-rays  of  the  jet  at  early  times  were  then 
taken  in  subsequent  test  firings  which  confirmed  the  theoretical 
predictions.  Figure  21  shows  the  jet  obtained  from  one  of  the 
experimental  test  firings.  The  equatorial  protion  of  the  liner  is 
seen  in  the  process  of  pinching  off  the  jet. 
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V.  COMPARISON  OP  DIFPBRHNT  DBSIGN  CONFIGURATIONS 

The  above  numerical  study  dealt  exclusively  with  the  hemispheri- 
cal liner  collapse  under  implosion  and  the  jet  obtained  from  the 
shaped  charge  configuration  shown  in  Figure  22Ca).  In  this  section 
we  shall  compare  briefly  the  collapse  and  jet  properties  of  the 
charges  shown  in  Figures  22(b)  through  22(f)  with  the  above  studied 
charge  22(a). 

The  charge  shown  in  Figure  22(b)  is  cylindrical  in  shape,  contains 
the  same  mass  of  explosive  as  the  previous  case,  but  is  initiated 
at  its  end  by  a plane  wave.  Figure  23  shows  the  computational  setup, 
while  Figure  24  shows  the  plane  detonation  wave  at  t - 2.2ms  shortly 
before  striking  the  liner.  Figures  25  and  26  show  the  early  stages 
of  the  liner  collapse  in  this  case.  Figures  27  and  28  show  the  early 
stages  of  jet  formation  for  this  case,  while  Figures  29  and  30  show 
the  late  stage  jet  configurations. 

It  is  apparent  from  examining  Figures  23  through  30  and  the  rest 
of  the  computational  plots  (not  shown)  that  the  liner  turns  inside 
out  In  this  case.  The  resulting  jet  from  this  shaped  charge  design, 
which  is  defined  to  be  that  part  of  the  collapsed  liner  with 
velocities  ^ 2mm/us,  is  seen  to  be  a short  thick  jet.  Some  instability 
is  seen  to  occur  near  the  base  of  the  slug.  The  jet  takes  a longer 
time  to  form,  its  mass  is  about  the  same  as  the  previous  case  22(a) 
(11.6  xl0“2kg),  but  is  considerably  slower  moving,  the  jet  tip  velocity 
being  only  (3.6mm/MS),  The  total  energy  communicated  to  the  liner  is 
(8.9  X lO^J)  which  is  about  52%  of  the  energy  communicated  to  the 
liner  in  tho  previous  case  22(a).  Figure  31,  views  (a>  through  (c) , 
show  the  jet  obtained  experimentally  from  a copper  hemisphere  whose 
radius  is  1.9  cms  with  a cylindrical  charge  which  is  point  initiated 
at  its  end  along  the  axis  of  symmetry.  Figure  31(d)  shows  our 
computational  prediction  of  the  liner  configuration  at  t = 28.7ms 
after  initiation  for  the  charge  shown  in  Figure  22(b)  with  a plane 
wave  as  described  above.  The  similarity  between  the  experimental  and 
computational  results  is  encouraging  and  provides  some  validation  of 
our  computational  model. 

The  charge  shown  in  22(c)  is  similar  to  22(a)  except  that  some  of 
the  explosive  in  the  equatorial  region  was  shaved  off  in  the  hope  that 
the  equatorial  portion  of  the  liner  will  take  a longer  time  before 
it  interrupts  the  jet.  The  interruption  of  the  jet  was  delayed 
about  6ms,  but  it  was  not  possible  to  eliminate  it.  The  jet  obtained 
in  this  case  was  similar  to  that  obtained  from  22(a),  Table  1 contains 
the  characteristics  of  the  jets  obtained  from  charges  22(a)  through 
22(f),  while  Figures  32  through  35  show  some  profiles  of  the  jets 
and  liners  obtained  from  charges  22(c)  through  22(f)  respectively  at 
the  times  indicated  on  the  Figures. 
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Table  1.  Summary  of  the  Computational  Results. 
CHARGES  GEOMETRY 


FIGURE  22 

(a) 

(b) 

Icl 

Cd)* 

(a) 

CO 

Jet  Mass  (Kg) 

xlO^ 

11.7 

11.6 

11.8 

9.6 

12.6 

18.3 

Jet  Tip  Velocity 
(rom/ps) 

6.4 

3.6 

5.0 

6.1 

6.8 

5.7 

Kinetic  Energy 
of  Jet  (J) 

xlO'*^ 

6.36 

3.72 

5.81 

5.74 

8.35 

9.83 

Total  Liner 

Energy  (J) 

xio"^ 

1.68 

0.89 

1.63 

1.17 

1.94 

2.90 

The  charge  shown  in  Figure  22(d)  has  an  iron  flange  6.3S  mm 
(0.25  in)  thick  to  confine  the  expansion  of  the  detonation  products 
at  the  equatorial  surface  of  the  charge  and  delay  the  rarefaction 
emanating  from  that  surface.  The  results  of  that  computation  are 
summarized  in  column  4 of  Table  1,  however  those  numbers  listed  are 
believed  to  be  in  error  because  parts  of  the  jet  flowed  out  of  the 
computational  mesh  at  the  top  transraittive  boundary  of  the  computa- 
tional mesh  before  stablizing  into  a final  state.  Normally  such  a 
problem  is  avoided  through  constant  monitoring  of  the  calculations 
and  a timely  rezoning  of  the  computational  mesh.  The  numerical  results 
given  in  column  (d)  should  therefore  be  considered  as  lower  bounds 
for  this  case.  In  the  numerical  simulation  of  the  performance  of 
this  configuration,  the  detonation  gases  were  not  allowed  to  expand 
between  the  liner  and  the  flange  as  it  is  the  case  in  practice.  It 
was  observed  in  the  computational  model  that  the  flange  prevents  the 
jet  disruption  by  the  liner's  equator.  Figure  33  shows  the  flow 
field  configuration  at  t = 13.49ms. 

The  charge  shown  in  22(e)  is  similar  to  22(a)  but  has  a confining 
case  on  the  outer  surface  of  the  charge.  The  jet  obtained  in  this 
case  has  31%  more  kinetic  energy  than  the  unconfined  charge  22(a). 

The  jet  has  a larger  mass  and  a faster  jet  tip  velocity.  A larger 
portion  of  the  explosive  energy  is  communicated  to  the  liner  in  this 
case.  The  results  of  this  computation  are  summarized  in  column  S 
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Figure  25.  Early  Stage  of  Liner  Collapse,  Charge  22(b) 
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of  Table  1.  Unfortunately  the  equatorial  section  of  the  liner  Inter- 
rupts the  jet  at  an  early  stage  of  Its  fomatlon.  Figure  34  shows 
the  liner  configuration  at  t « 12.8Svis. 

The  charge  shown  In  22(f)  has  a 3.81  mm  (0.150  In)  thick  liner 
which  Is  double  the  thickness  of  the  liners  used  In  the  previous 
cases.  The  charge  thickness  was  also  Increased  to  24.13  mm  (0.95  in). 
The  Initiation  points  were  distributed  over  a 120”  sector  of  the  charge 
surface  Instead  of  the  180”  hemispherical  surface.  The  calculations 
revealed  that  the  equatorial  portion  of  the  liner  interrupts  the  jet 
about  ISps  after  initiation.  The  jet  has  more  mass,  but  the  jet 
velocity  is  slightly  lower  than  case  (a) . A summary  of  the  results 
in  this  case  is  given  in  column  6 of  Table  1,  and  Figure  35  shows 
the  flow  field  configuration  in  this  case  at  t » 13.02ps. 


VI.  CONCLUSION 

The  above  study  explained  the  collapse  and  jet  formation  processes 
of  a hemispherical  liner  with  a hemispherical  charge  which  is  surface 
imploded.  It  is  clear  from  the  preceding  study  that  the  collapse  and 
jet  formation  processes  are  fundamentally  different  from  those  of 
a cone.  The  various  hemispherical  liner  elements  converge  towards  a 
single  point  as  they  collapse.  The  jet  forms  due  to  the  compression 
of  liner  material  which  squeezes  out  the  jet  materials  as  in  an  extrusion 
process.  Conical  liners  collapse  along  a line.  The  study  has  shown 
that  the  jet  originates  from  the  inner  surface  layer  of  liner  material. 

It  has  been  shown  that  such  jets  have  no  inverse  velocity  gradients. 

The  above  study  revealed  the  cause  of  the  jet  pinchoff  as  being  due 
to  the  equatorial  portion  of  the  liner  striking  the  jet.  The  design 
modifications  discussed  in  Figures  22(c)  and  22(d)  alleviate  this 
problem. 

The  second  part  of  the  study  dealt  with  the  effects  of  variations 
of  the  charge  geometry  and  initiation  mode  as  in  charge  22(b)  and 
dealt  with  the  effects  of  confinement,  and  variations  of  the  charge 
and  liner  thicknesses.  The  study  of  charge  22(b)  revealed  that  the 
hemispherical  liner  turns  inside  out  in  this  case  and  forms  a short, 
thick,  and  slow  moving  jet  which  might ^be  suitable  fox  use  against  a 
certain  class  of  targets. 

In  conclusion  the  above  computational  study  has  provided  a wealth 
of  quantitative  values  of  various  physical  parameters  most  of  which 
agree  quite  well  with  experimental  measurements.  This  agreement  in 
addition  to  the  qualitative  agreement  of  our  computations  with 
experimental  results  provides  validation  of  our  computational  model 
and  reveals  the  benefits  of  such  studios  as  a design  tool.  In 
practice  many  expeximontal  design  parameters  were  based  on  similar 
hydrodynamic  computations.  Such  computational  studies  can  often  be 
made  quickly,  economically  and  can  result  in  considerable  financial 
and  time  savings. 
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